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Abstract. We extend a result of Stolz and Weidmann on the approximation 
of isolated eigenvalues of singular Sturm-Liouville and Dirac operators by the 
eigenvalues of regular operators. 



1. Introduction 

The approximation of isolated eigenvalues of singular ordinary differential oper- 
ators by the eigenvalues of regular operators is an important and well studied topic 
since the latter ones can be computed numerically with arbitrary precision. See the 
recent monograph by Zettl [5] or in particular the recent survey [7J by Weidmann. 

While the case of eigenvalues below the essential spectrum is well understood, 
the case of eigenvalues in essential spectral gaps was only recently solved by Stolz 
and Weidmann in [3] (see also [1]). 

Let r be either a Sturm-Liouville expression 

, 1 / d d 

(1) T=- (-—p— + 

r V dx dx 



or a Dirac system 

1 /. d 

(2) T=- R72 



dx 

where 02 = f ? is the Pauli matrix, on the interval / = (a, b). 

As usual, we will assume the coefficients p^ 1 , q, r are real- valued locally integrable 
functions with p, r > in the Sturm-Liouville case and q, r are real, symmetric 2x2 
matrices with r > in the Dirac case. 

An endpoint a or & is called regular if it is finite and the coefficients are integrable 
near this endpoint. If both cndpoints are regular, we will call r regular. 

Let 33 (t) be the maximal domain of definition of r. Then r is sclf-adjoint on 
23 (t) if it is limit point (l.p.) near both a and b. Otherwise we will impose an 
additional boundary condition at every endpoint where r is limit circle (I.e.). In 
this way we obtain a self-adjoint operator H associated with r. 
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A solution ip a (z, x) (ipb(z,x)) of rip — zip which is square integrable near a (b) 
and satisfies the boundary condition at a (6) (if any) is called a Weyl solution. Such 
a solution is unique up to a constant and it exists at least for z 6 C\cr ess (H). 

Our aim is to approximate H by regular operators H n obtained by restricting r 
to a finite interval (a n ,b n ) C (a, 6) (O Chap. 14]). The case a n = a or b n = b is 
allowed. 

Fix functions u, v G £)("?")• Pick a n J, a, 6„ f b. Define H n 

(3) #„ : D(H n ) -> L 2 ((a„,6„);rdx) , 

/ -> r/ 

where 

(4) D(ff n ) = {/ e S(rn))|W „(tx,/) = W K (v,f) = 0}. 
Then Stolz and Weidmann prove the following 

Theorem 1 ([3]). Define H n as above with u — ip a {X a ) and v — ipb(Xb) with 
X a ,Xb € [Ao,Ai] (in particular, assume that the corresponding Weyl solutions exist). 
Let Pfi(H) be the spectral projection of H corresponding to the Borel set fl C R. 

// dimRanP(^ ^ t \(H) < oo, then the eigenvalues of H in (Ao,Ai) are exactly 
the limits of eigenvalues of H n which lie in (Ao,Ai). The corresponding (one- 
dimensional) eigenprojections converge in norm. 

If dimRan_P(^ ^(H) = oo, then the eigenvalues of H n accumulate in (Ao,Ai) 
as n — > oo. 

In fact, in [3] this result is only proven in the case [Ao, Ai] na ess (H) = and Ao, Ai 
are not eigenvalues of H. Hence we will provide a proof of this slightly generalized 
version below. 

As pointed out in [3], this result has of course one practical drawback: The 
Weyl solutions used to generate the boundary conditions of H n will not be known 
explicitly in general. To evade this obstacle they show that their result still holds 
if the Weyl solution of a nearby operator is chosen instead. 

The main purpose of this note is to propose an alternate way of approximating 
H which only involves one Weyl solution at one endpoint. More precisely we show 
that if r is l.p. at one end point, the Weyl solution of the other endpoint can also 
be chosen instead: 

Theorem 2. Suppose r is l.p. at b. Define H n as above with u = ^a(A a ), 
A a S [Ao,Ai], and v = ipa(^o) or v — V'a(Ai) (in particular, we assume that the 
corresponding Weyl solutions exist). 

If dim Ran ^(H) < oo, then the eigenvalues of H in (Ao,Ai) are exactly 
the limits of eigenvalues of H n which lie in (Ao,Ai). The corresponding (one- 
dimensional) eigenprojections converge in norm. 

If dimRanP(^ ^(H) — oo, then the eigenvalues of H n accumulate in (Ao,Ai) 
as n — > oo. 

In particular, if one endpoint is regular, a solution satisfying the boundary condition 
at this endpoint can be taken in this case. Clearly the result cannot hold if r is I.e. 
at b, since our assumptions contain no information on the boundary condition of 
H at b in this case. 

Remark 3. (i). As shown in [3], if t is I.e. at a, then u can be choosen to be any 
function in D(t) generating the boundary condition of H at a. 
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(ii). The same result (with the same proof) holds for Jacobi operators (see ^\). 
2. Approximation by regular operators 

We begin by recalling that H n converges to H strongly ([6j). Strictly speaking 
this statement makes no sense since H n and H live in different Hilbert spaces. This 
can be easily fixed by using at © H n © al on L 2 ((a, b);rdx) — L 2 ((a, b n );rdx) © 
L 2 ((a n , b n ); r dx) ©Z/ 2 ((a„, 6); r dx) where a is a fixed real constant outside [Ao, Ai]. 
Alternatively, one can also use generalized strong convergence as introduced in [3J . 

Lemma 4. Suppose that either H is limit point at a or that u — i/>_(Ao) for some 
Ao and similarly, that either H is limit point atb or v = Tp + (Xi) for some Ai. Then 
H m converges to H in strong resolvent sense as m — > oo. 

In addition, we need the following abstract result. In this respect we remark that 
for a self-adjoint projector P we have 

(5) dimRan(P) = tr(P) = ||P||i, 

where ||.||i denotes the trace class norm. If P is not finite rank, then it is of course 
not trace class and all three numbers are equal oo. 

Lemma 5. Let A n , A be self- adjoint operators such that A n — > A in strong resol- 
vent sense. Then 

(6) tr(P (Ao>Al) (A)) < liminf tv(P (XoM) (A n )). 
If in addition. tr(P (AoiAl) (A„)) < tr(P (Ao:Al) (A)), then 

(7) lim tr(P (Ao , Al) (A„)) = tr(P (Ao , Al) (A)) 

and z/tr(P(x ,Ai)(-^)) < oo we even have 

(8) hm ||i > (A Oj A 1 )(^»)-- p {Ao,A 1 )(^)lli = 0- 

n — *oo 



Proof. The first part is just Lemma 5.2 from pQ. This also implies the second if 
tr(P( A()jAl )(^4)) = oo. Otherwise, if tr(P( AoiAl )(A)) < oo, we have 

limsuptr(P (AoiAl) (AO) < tr(P (Ao , Al )(A)) 

and the first claim follows. The second is then a consequence of Grumm's theorem 
(0 Thm. 2.19]). □ 

Now it remains to show that this result is applicable in our situation. 

Lemma 6. Let H n be defined as in (0) with 

(i) u = ipa(K)> v = ^&(Ab) with A a , A 6 £ [A , Ai] or 

(ii) u = ip a (X a ) with A Q € [A , Ai] and v = ^ a (A&) with X b E {A , Ai}. 
Then, 

(9) tr(P (AojAl) (P„))<tr(P (Ao , Al) (P)). 

Proof. Abbreviate P = P (Ao , Al) (P), P„ = /' \ , A .(//„). 

(i). Since this part is identical to the proof in [3], we just give an outline. Let 
ipi , . . . , ipk G Ran P n be the normalized eigenfunctions of H n , construct 

Ta,j , x <C a n , 
i'j( x )={ ^jO* 1 )) a n <x<b n , 
j b jv(x), x > b n , 
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where J a ,j, 7fc,j are chosen such that tpj £ 2(t). A computation now shows that 

for any ip in the linear span of the V'j's, which yields the first result. 

(ii) By considering two steps from (a n ,b n ) to (a, b n ) and from (a,b n ) to (a, b), 
we see that the first step is covered by (i) and hence it is no restriction to assume 
a n = a. Now proceed as in the previous case but use 



( x ) = | ~ iMx), X < b n , 



^ 1 0, X>b n , 

where jj are chosen such that ipj 6 5) (r). Now let ^ = ■ Cjtpj be in the linear 
span of the ipj's. Then, since v is also an eigenvector of H n and hence orthogonal 
to the "0j' s j we have 

ii trj ^i + > , ,|2 n /2Aj - Ai — A ~ , . 2A 6 — Ai — Aq , . 

II i H 7, — ml = II z^ c n — o ^'C*) 9 t^O*) 



-/^l c il V ^ ; ~ r|n v 2 y l|V|l (o.6n) 

3 

where 7 = Y]j Cjjj- Hence the second result follows. □ 
Theorem [T] and Theorem [5] now follow by combining the last two lemmas. 
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